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1. INTRODUCTION 

Portfolio optimization is concerned with making a balance between risk and return of the selected 
securities in the security market. Thus, the different composition of securities in the portfolio will give varied 
results. Due care must be taken in selecting right combination of securities to make the optimal portfolio. The 
very basic principle of investment is diversification where investors have to commit their fund into different 
securities/assets based on their respective return and risk. The diversification of the portfolio is not as simple 
as people think by committing funds into different classes of assets. Systematic diversification requires a 
number of inputs for different securities like their expected return, standard deviation of the return, variance 
and covariance of the return, coefficient of correlation between the returns of different assets. The 
diversification will lead to optimize risk and return [1], [2]. 

In the field of portfolio, the selection of best combination of securities out of the available in the 
security market has always been a tedious job for the experts of the field. Henry Markowitz has developed a 
model for portfolio optimization by emphasizing on selection of best securities in the portfolio. He 
emphasized on the computation of risk and the return of different securities and based on that optimal 
number of securities was selected to make the portfolio to its optimum level. In the real world, sometimes 
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returns of different securities are not symmetrically distributed but the model considers that the returns are 
normally distributed [3]-[5]. There have been many models suggested for the same purpose, but the model 
discussed by Markowitz has been the core model and inspires many others in the domain [1]. However, the 
Markowitz model for portfolio optimization may be suitable if number of variables is few. But in case of 
large number of variables or constraints, this model may not give authentic and reliable result. Due to this 
limitation of the model suggested by Markowitz, some other techniques have been developed to give better 
solution by mixing the technique with quadratic programming. But attaining the optimum portfolio will be 
more difficult if the additional constraints are considered like boundary constraints, cardinality constraints 
which are in the form of non-linear mixed integer programming problems. 

Solving such kinds of problems are extremely difficult than original problems. The available 
solutions are not enough to handle such types of problems. In such cases, the swarm intelligence (SI) and 
evolutionary computation (EC) approaches are being used in construct optimal portfolio by predicting the 
global optimum. Chang ef al. [6] have been reported that genetic algorithm (GA) based solution for portfolio 
optimization problems suitable to approximate the unconstrained competent frontier. Kyong et al. [7] and Lin 
and Liu [8] presented the GA to optimize for the index fund management and transaction costs, respectively. 
Cheng et al. [9] proposed a novel approach for portfolio selection problems-based GA to measures the three 
different risk such as semi-variance, variance with skewness and mean absolute deviation. Jalota and Thakur 
[10] have designed to solve the portfolio selection problems by GA for handling the cardinality constraint, 
lower/upper bound and budget constraints. Li [11] presents a novel approach based on GA to solve the 
problems of investment and income for enterprises or individuals by mixing of operations research and 
finance. Chang and Hsu [12] proposed a particle swarm optimization (PSO) algorithm to select the top five 
portfolios of the stocks from the equity fund to optimize the ratio return and return rates. Jiang et al. [13] 
proposed PSO algorithm based on diffusion repulsion to portfolio selection problem to keep the faster 
convergence rate. Zhu ef al. [14] proposed a novel approach using PSO for portfolio optimization problems 
to test the unrestricted and restricted risky portfolio investment. Reid and Malan [15] applied PSO algorithm 
to handle the two new constraints likes preserving feasibility and portfolio repair method. Deng and Lin [16] 
applied ant colony optimization (ACO) for mean-variance portfolio optimization model to solve the 
cardinality constraints to effective low-risk investment. Bacterial foraging optimization (BFO) [17] is another 
powerful approach for optimization problems. Niu et al. [18] proposed a BFO algorithm for portfolio 
selection problem to optimize liquidity risk by introducing the endogenous and exogenous liquidity risk. 
Kalayci et al. [19] proposed an effective solution approach based on an artificial bee colony (ABC) algorithm 
with infeasibility toleration procedures and feasibility enforcement for solving cardinality constrained portfolio 
selection model with the aim to optimize the return of investment. Chen et al. [20] reported an improved version 
of ABC algorithm for portfolio optimization to focus on balance the trade-off between return and risk. Mazumdar 
et el. [21] proposed a novel approach for portfolio and unsystematic risk selection problem using grey wolf 
optimization to minimize the risk contributor that improve the diversification ratio. Shahid et el. [22] 
proposed a gradient based optimizer for unconstrained portfolio selection model. In another work [23], an 
invasive weed optimization has been applied for risk budgeted portfolio selection model optimizing Sharpe 
ratio. 

Some hybrid approaches are also reported to combine meta-heuristics with exact algorithms or other 
meta-heuristics to deal with complicated portfolio selection problem. Maringer and Kellere [24] used a 
hybrid local search algorithm that combines a meta-heuristic namely simulated annealing (SA) and 
evolutionary algorithms (EA) for cardinality constrained portfolios. Tuba and Bacanin [25] proposed a hybrid 
approach combination with ABC and e firefly algorithm (FA) to optimize the mean variance return, 
Euclidean distance and retune error. Qin et al. [26] proposed a novel hybrid algorithm based PSO and ABC 
for conditional value at risk of portfolio optimization problem to optimize the mean and standard deviation. 

In this paper, authors proposed a portfolio selection model to maximize the Sharpe ratio by using 
stochastic fractal search (SFS) based evolutionary optimization approach. This approach is derived from 
natural growth process which is mathematically modeled by fractal theory to explore the solution space with 
a number of constraints. The major contributions of the work are listed as: 

— To apply SFS based evolutionary optimization approach maximizing Sharpe ratio of the portfolio 
constructed. 

— Performance comparison has been done with state of art population-based models (GA, SA) [6] from 
the domain. 

— Experimental analysis has been conducted by using real datasets of the Bombay stock exchange (BSE). 

The rest of the paper is organized as: section 2 presents the mean-variance model for Sharpe ratio 
with constraints handling procedures. Further it describes the SFS approach with algorithmic template. 
Section 3 reports the experimental results and their interpretations. Section 4 draws conclusion of the paper. 
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2. RESEARCH METHOD 

In this section, the problem statement with a mathematical model of the portfolio selection model 
has been presented. Proposed evolutionary algorithm-based solution approach has been discussed in detail. 
Further, an algorithmic template has been also given for better understanding of the proposed solution 
approach. 


2.1. Problem formulation 

In investment market, a portfolio (P) is designed with K assets from asset set, A = { a;,Q2,..ax} 
with weight set, W = { 1, W2, ... Wx}. Their respective expected returns are represented by using a return 
set, such as, R = { Ty, 12, ++ Tx}. For the scenario, total portfolio risk and return can be expressed by using (1) 
and (2): 


Prisk = i Li wu; *w; * CoV(i,j) (1) 


Preturn = par yew (2) 


where w;, and w; are the weights of a; and aj, respectively. CoV (i,j) is the covariance matrix generated by 
return values of the assets in specified duration. 

In portfolio design, investors target to maximize Sharpe ratio(SRp) of the portfolio under 
consideration. Risk free return is assumed here as zero as we are taken equity-based assets. Now, the Sharpe 
ratio of the designed portfolio is estimated and can be written as (3). 


Max (SRp) = max (Feu) (3) 


Prisk 


Subject to the constraints 
K = 
a) Met WwW; =1 


c) a Su, <b; 

Here, (a) represents budget constraint. Further, (b) constraint restricts the short sell. Lastly, (c) 
constraint imposes lower and upper bounds for assets weights. Constraints listed above from (a) to (c) are 
repaired by using the constraint handling procedures given as follows: 

Moreover, aS we see here, the above constraints from (a) to (c) are linear with convex feasible 
region. Therefore, a generalized penalty method used for constrain handling is also presented in this section. 
Now, consider, the search space of the decision variables is represented by X, then the penalty function M 
may be as (4). 


_ (G+ Vi gi), ifx EX 
i leas if xX @) 


where, G and g are the solution and the constraints value of x in the feasible space. 


2.2. Stochastic fractal search (SFS) 

In this section, solution approach portfolio selection using SFS has been presented to make effort for 
optimal weights to maximize Sharpe ratio. This algorithm has been proposed by Salimi [27] which is an 
evolutionary population-based method with two phases, namely, the diffusion and updating. It is derived 
from the natural growth modeling by using concepts of fractal theory. In diffusion, diffused points create new 
points by moving around neighboring positions, also avoids premature convergence by avoiding local 
optima. The diffusion phase is presented by applying any of two Gaussian walks (G;): 


G, = Gaussian(t1p,,.,.0) + (5 * Poest — 8’ * Pi) G) 
G2 = Gaussian(Up, 0) (8) 


Here, P; is the location of i-th point, best point (Py.e<,) represents the best position of iteration, and 6, 6’ € 
[0, 1]. The parameters involved in the Gaussian function (up,,., and Up) are same as |Pyest| and |P;|. And, 


er F : 1 
the standard deviation (0) is estimated as o = —_ * (P; — Prest)]- 
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Next, the updating phase alters every point’s positioning as per the most suitable position of a point 
in the population. Two update procedures are as: 


PiG) = P.G) — € * (RW — YiG)) (9) 
Pi —6«(PL-P if 6 <0.5 
Pr’ = 7 a ( g best) (10) 
P/+8+(Pi-Pi) if 6’ >05 


In (9), P. and P; are the two randomly generated distinct points obtained from diffusion phase. And, in (10), 
P; is the new point obtained from (9). And, 6 & 6 € [0,1], are random numbers. Finally, the proposed SFS 
algorithm to get the optimal weights for considered portfolio has been given as follows: 


Algorithm : SFS () 
Begin 
1. Initialize N points population 


2. while (i < largest generation) do // i: number of iterations 
2% for each P; do // PP; : i Points from N points 
4. Diffusion () 
Die { 
6. G=" Gax // maximum number of diffusions 
Ts for j = 1: q do 
8. if (G, is selected) 
9. generate a new point by using (7) 
0. else if (G, is selected) 
ale generate a different point by (8) 
Qi end if 
Se end for 
4. } 
on end for 
6. Update-I () 
7. { 
8. Rank the points as per fitness values 
Oi for each P; do 
20 for each j —P; do 
21 if (rand [0, 1] 2 Pa;) 
22. update jth component of P; by (9) 
23% else if 
24. do nothing 
25:. end if 
2-6. end for 
27. end for 
28%. } 
29. Update-II () 
30 { 
31. Rank all points from update-I on fitness values 
32%. for each new P; do 
33% if (rand [0, 1] 2 Pai) 
34. update the point by (10) 
35% end if 
36. end for 
B's } 
38. end while 


End 


3. RESULTS AND DISCUSSION 

In this section, a performance evaluation has been planned to evaluate the proposed model with 
comparative experimental analysis. This analysis has been done on an Intel (R) Core (TM) i7 CPU 3.20 GHz 
processor with 16 GB of RAM using MATLAB. Some well-known meta-heuristics namely genetic algorithm 
(GA), simulated annealing (SA), are also considered for comparative analysis. In the experimental analysis, 
datasets (BSE 30, BSE 100, BSE 200, and BSE 500) are extracted from S&P BSE Sensex of Indian stock 
exchange of monthly holding period returns from 1* April 2010 to 31 March 2020. The parameter setting is 
more challenging issue in meta-heuristic approaches. In SFS, Size of initial population (n) and Maximum 
number of iteration (Itermax) control the rate of convergence. Iter a, can be set as per the solution 
optimality requirement, a larger value of Iter,,,, makes the method able to achieve a better result. 
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In this work, the best combination of parameters for the algorithm is used determined for 
experiments. Thus, the proposed single-objective model optimizing Sharpe ratio is solved with best set of 
parameters on all considered datasets considering other parameters fixed. The parameters for proposed and 
other algorithms for comparative analysis are listed in Table 1. 


Table 1. Control parameters for GA, SA and proposed SFS algorithms 


Algorithms Parameters Specifications 
Common parameters n= 100, Iter,,,,= 100 
SFS Random walk =1, d = 25, Pa; = 1/2 
GA Simulated crossover probability = 0.7 
Polynomial mutation probability = 0.3 
SA Neighbors =5, Mutation probability = 0.5 


Here, Figures 1(a) to (d) demonstrates a convergence comparison of GA, SA and SFS algorithms 
over objective fitness function. The bench mark datasets from Bombay Stock Exchange namely BSE 30, 
BSE, 100, BSE 200 and BSE 500 has been used for the study. For performance comparison, simulation 
experiments of each algorithm Viz. GA, SA, and SFS has been taken in the study. 


BSE 30 BSE 100 
0.4 ; 


0.3 


Best fitness function values 
Best fitness function values 


0.25 fi ‘ 
0 20 40 60 80 100 20 40 60 80 100 
Iteration Iteration 
(a) (b) 


BSE 200 BSE 500 


Best fitness function values 
Best fitness function values 


0 20 40 60 80 100 
Iteration Iteration 
(c) (d) 


Figure 1. Convergence curve of GA, SA, and SFS on; (a) BSE 30, (b) BSE 100, 
(c) BSE 200, and (d) BSE 500 


We have performed 20 different runs to obtain the best optimal value of weight vector maximizing 
the Sharpe ratio for all data sets in the formulated optimization model. And the Max (best), Min (worst) and 
average values of Sharpe ratio of obtained solutions by various models are reported in Table 2. The best 
values among the considered models are shown in bold in the tables. The results of GA, SA, SFS algorithms 
of various runs for Sharpe ratio of all datasets are presented in Box Plots in Figures 2(a) to (d) for better 
presentation and graphical self-interpretation. 
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Table 2. Sharpe ratio 


GA SA SFS 
K=30 Max 0.390334 0.396953 0.400857 
Min 0.365241 0.394756 0.400597 
Avg. 0.372709 0.395917 0.400829 
K=100 Max 0.471465 0.439110 0.486310 
Min 0.399370 0.425896 0.471068 
Avg. 0.440248 0.430771 0.482106 
K=200 Max 0.556202 0.477881 0.604232 
Min 0.450313 0.435389 0.544238 
Avg. 0.522603 0.456614 0.575661 
K=500 Max 0.608987 0.378147 0.636084 
Min 0.462989 0.298987 0.542441 
Avg. 0.535208 0.351277 0.590286 
BSE30 BSE100 
0.4 fK = 
== 0.48 F = 
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Figure 2. Box plots of all four algorithms on (a) BSE 30, (b) BSE 100, (c) BSE 200, and (d) BSE 500 


SFS algorithm performs better in account of both the convergence rate and fitness value than GA 
and SA as shown in Figures 1(a) to (d). Further, SFS performs remarkably good compared to others and 
delivers higher values of Sharpe ratio (on all measures i.e, Max., Min., and Avg.) for all considered 
benchmark datasets such as BSE 30, BSE, 100, BSE 200 and BSE 500. The performance order on objective 
parameter of remaining algorithms is GA and SA. Thus, we can argue that the proposed approach greatly 
contributes robust portfolio optimization with satisfied desired constraints. 


4. CONCLUSION 

The portfolio selection problem is one of the core problems in investment management which drawn 
the attention of investors in the recent decades. Due to constraints need to be managed in portfolio 
construction, the conventional techniques are not good enough in giving solutions. Therefore, recent 
optimization methods are used to find optimum solution in complex scenario. In this work, a portfolio 
selection model based on an evolutionary algorithm namely SFS has been proposed. The natural growth has 
been mathematically modeled to explore the search space for optimum solution. For performance evaluation, 
an experimental analysis has been conducted to determine the effectiveness of the proposed model by doing 
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performance comparison with state of art models from the domain such as GA and SA. The real datasets of 
the S&P BSE Sensex of Indian stock exchange have been taken for performance evaluation. Study shows the 
superior performance of SFS on objective parameter among its peers for all data sets in the study. 
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